Coulomb blockade in graphene nanoribbons 
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We propose that recent transport experiments revealing the existence of an energy gap in graphene 
nanoribbons may be understood in terms of Coulomb blockade. Electron interactions play a decisive 
role at the quantum dots which form due to the presence of necks arising from the roughness of the 
graphene edge. With the average transmission as the only fitting parameter, our theory shows good 
agreement with the experimental data. 
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Graphene, a two-dimensional allotrope of carbon on a 
honeycomb lattice, was isolated a few years ago [1\ cre- 
ating a great excitement in the physics community due 
to its close connections to high-energy particle physics 
[3, 0] and its tantahzing possible technological applica- 
tions 0, 0| • It is now experimentally established that a 
great deal of the properties of graphene ^ can be de- 
scribed in terms of non-interacting (or weakly interact- 
ing) linearly dispersing Dirac quasi-particles [3, [11 . The 
only accepted exception maybe when graphene is subject 
to strong magnetic fields in the quantum Hall regime (sf , 
when the electronic kinetic energy is quenched by the 
appearance of Landau levels, and the fourfold degener- 
acy of the Landau levels is split by electron-electron in- 
teractions [13, Hi- Nevertheless, in the absence of an 



applied magnetic field, because of the vanishing of the 
density of states for Dirac fermions in two dimensions 
[l^ . the electrons in graphene interact through strong, 
essentially unscreened , long-range Coulomb interactions 
I3L 14, 3, 3]- The fact that Coulomb interactions do 
not show up in bulk experiments remains a puzzle in the 
physics of graphene. 

Recent experiments on the electron transport proper- 
ties of lithographically patterned graphene nanoribbons 
have shown the existence of an energy gap near the 
charge neutrality point [13] ■ The size of the gap Eg is 
inferred from the nonlinear conductance at low temper- 
atures and is found to decrease with the ribbon width 
W following the approximate law Eg = a/{W — W*). 
This result seems to correlate with a conductance behav- 
ior G = /3{W - Wo), since M/* « Vl^o - 16 nm for the 
same sample at temperature T = 1.6 K. In the absence 
of interaction effects, the energy gaps between subbands 
in a graphene ribbon should scale inversely with the rib- 
bon width. Eg w hvp/W, where vp ~ 10^ m/s « 0.66 eV 
X nm X is the Fermi velocity in graphene [il]. Nev- 
ertheless, this estimate leads to gaps which are smaller 
than those observed experimentally (note that, for widths 
of about 20 nm the experimental gaps are larger than 
0.1-0.2 eV). This result has led to the suggestion that 
the effective transport width is reduced with respect to 
the nominal width W by an amount W* « Wq due to 



the existence of structural disorder at the edges or to a 
systematic inaccuracy in the determination of the geo- 
metrical width caused by over-etching beneath the etch 
mask [l3|- On the other hand, it has been shown that 
graphene quantum dots as large as 25 /im (at low temper- 
atures), and as small as 40 nm (at room temperature), 
show Coulomb blockade effects [6| indicating that elec- 
tron interactions become stronger as the dimensions of 
graphene sheets are reduced. 

In the present work we argue that the main results of 
Ref. [17[ can be naturally explained as due to Coulomb 



blockade effects originated by the roughness at the edges 
of graphene nanoribbons. This roughness occurs natu- 
rally in graphite samples, as has been seen in scanning 
tunneling microscopy [l8| . leading to the localization of 
char ge a t the edges and the formation of electronic pud- 
dles [19| . As depicted in Fig. [U disorder at the edges of 
a graphene ribbon also leads to the formation of "necks" , 
causing an abrupt reduction in the number of conducting 
channels and thus to a large increase in the impedance 
along the graphene sheet. This results in the electric 
isolation of nanoscale size regions, or "dots" , where the 
electrons become temporally confined. Within this pic- 
ture. Coulomb blockade [H [HI, [H, [H, [Hi results from 
the electron transport from dot to dot through graphene 
necks. 

In our analysis, we will follow the standard theory of 
Coulomb blockade effects, neglecting the electronic level 
spacing within the grains. Note that graphene states 
delocalized throughout a region of linear size W show 
level spacings of order hvp/W, which is the same scaling 
behavior followed by the charging energy of a grain of 
size W, e^/W. It seems likely, however, that the rough 
edges of the samples studied in [l^, as well as the in- 
ternal lattice defects, can lead to a variety of partly or 
fully localized states at energies close to the Fermi level, 
reducing the electronic level spacing [1, [25, 2^, 13]. 

In the presence of Coulomb blockade characterized by 
a renormalized charging energy E* (which here we as- 
similate to Eg), the conductance between neighboring 
metallic dots is renormalized to lower values as the en- 
ergy or temperature scale is reduced [H, [2^ [s^l ■ This 
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FIG. 1: Illustration of a graphene ribbon with a disordered 
edge leading to the formation of necks and dots along the 
ribbon. Coulomb blockade takes place when the charge moves 
from dot to dot. 



renormalization begins to be appreciable at temperatures 
T < Ec/ks, where Ec is the charging energy associated 
to the geometric capacitance of the dot. A perturbative 
analysis, valid for relatively high temperatures such that 
T ~ Ec/ks, leads to a conductance given by: 



(1) 



where c = 2 and Go — 'iGqNchT is the non-interacting 
conductance (we include spin and valley degeneracy), 
Gq = jh = 38.7 /iS being the quantum of conductance, 
A^ch = kpW ji: the number of transverse orbital channels, 
and r the average transmission per channel. Eq. ([1]) is 
valid when the transmission per channel is low, r 1, 
although the total conductance can take arbitrary values. 
A similar expression can be obtained in the limit r < 1, 
provided that none of the channels has perfect transmis- 
sion 3l|. In the latter case, the constant c in Eq. ([T]) 



takes the value c = S/tt . Within the relatively narrow 
range of widths and temperatures considered in Ref. [l7| , 
the logarithm in Eq. ([T]) can be taken as a constant of 
order unity, leading to the approximate expression: 



G~ (4GQfcWr)(W^- Wo) 



(2) 



where Wo ~ 'Kc/AkpT. 

The experimental dependence of G on W can be used 
to estimate the average transmission r. Figure 2 of 
Ref. [r^ indicates that, at room temperature, and for 
a gate voltage Vg — liiirac = —50 V, a change in con- 
ductance AG « 80 /iS takes place if the width changes 



by AW « 40 nm. This yields Wo « 40 nm. The corre- 
sponding hole density in this experiment is n = 3.6 x 10^^ 
cm^^, which implies that k^^ — Xjix^fn k, 1.7 nm. These 
experimental results are consistent with Eq. ^ if the 
average transmission per channel is r « 0.07. The to- 
tal number of channels, for a width « 40 nm, is 
4A^ch ~ AkpW/iT « 30. The average transmission prob- 
ability found here is consistent with tight-binding calcu- 
lations for wedge shaped graphene constrictions [32| (see 

also [allai). 

The analysis which at high temperatures leads to Eq. 
([J) also shows that, at low temperatures, the effective 
charging energy is renormalized by virtual jumps of the 
electrons across the junction, leading to: 



Eg w Ece-'^^''^'^ 



(3) 



where E^ is the charging energy for the completely iso- 
lated dot. It seems reasonable to assume that Ec is de- 
termined by the nominal ribbon width W , Ec ~ e^/W 
since that is the size of a typical puddle which is isolated 
from the rest of the ribbon through a contact (see Fig. [1]). 
This contact acts as the bottleneck which determines the 
conductance of the graphene nanoribbon. This fact does 
not preclude, however, the possibility of further structure 
in the I-V characteristic which could be induced by the 
presence of other dots. 

The capacitance of the grains is also modified by the 
presence of metallic leads and gates. The leads, and 
the regions of the quasi~one-dimensional ribbon at dis- 
tances greater than W from the island considered do 
not change appreciably the charging energy, as a one- 
dimensional charge distribution does not screen an elec- 
trostatic potential. The screening of a metallic gate 
at a distance d from the island can be analyzed, when 
d ^ W , by assuming that the charge in the island in- 
duces an image charge. The charging energy is changed 
to Ec ^ e^/W — e^/2d. In the following, we neglect the 
second term as ~ 20 — 100 nm and typical distances 
to the gate [l| are d ~ 300 nm. Charging effects should 
be strongly suppressed when the distance to the metal- 
lic gate is comparable, or smaller, than the width of the 
ribbon. 

Finally, we have identified the value of the Coulomb 
gap at low temperatures with the gap in the TV charac- 
teristics measured in Ref. 17[. Note that the argument 
in the exponent in Eq. ^ becomes 4/cf Wt/ttc ~ 0.8 
for W « 40 nm, which leads to an appreciable renor- 
malization of the geometrical charging energy by virtual 
charge fluctuations (see [IH] for similar effects in a dif- 
ferent granular system). We note that Ec ^ A meV for 
ly sa 40 nm. 

We can also write Eq. ([3]) as: 



Eg{W) 



W 



(4) 
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FIG. 2: Comparison between experimental data of Ref. [13] 
and the theoretical result (full line): = BWe'^^, ob- 

tained from Eq. Q with B = 0.001 (meV x nm)~^ and 
C = 0.023 nm"\ 



where Wq is the length scale used in Eq. Thus, as a 
function of W, E^^ takes low values for W ^ Wq while 



experiencing a sharp rise for W > Wq- We note that, 
E~^{Wo) « 5 X 10"'^ (meV)~-^, substantially smaller that 
the unit scale used in Fig. 3e of Ref. [l3| ■ One can see 
from Fig. [2] that Eq. ([4]) explains the experimental data 
of Ref. [17.] within its error bars. A comparison between 
the expression E~-^ = BWe^^ , obtained from Eq. (g]) 



and the data in 17| yields — 43 nm, which is in 

good agreement with our theoretical estimate for Wq. 
One also obtains B = 10^^ (meV x nm)^^, which agrees 
reasonably with — 6.95 x 10^^ (meV x nm)~^. The 
fact that is smaller than can be partly attributed 
to the screening effect of the gate. 

Finally, we note some trends which provide additional 
qualitative support to the Coulomb blockade picture. 
Figs. 3a, 3c, and 3d of Ref. [13] show the differential 
conductance as a function of the bias and gate voltage. 
Electron-hole symmetry explains the symmetric behav- 
ior around a value of Vg which must be identified with 
the neutrality point. The maximum vertical width of the 
dark (low differential conductance) zone must be identi- 
fied with the gap Eg{W), which clearly decreases with 
W, in qualitative agreement with Eq. ((4]). On the other 
hand, varying Vg is equivalent to changing k^- Eq. ^ 
shows that the variation of the gap and the differential 
conductance as a function of k-p must be faster for large 
values of W, in good agreement with the experimental 
results. 

The experiment of Ref. [13] also shows that the linear 
(Vf, 0) conductance depends weakly on the gate volt- 
age at the neutrality point, where kp — )■ 0. Our analy- 
sis predicts that the linear conductance depends on gate 
voltage through the product kpT. The well-known ex- 
istence of a minimum in the bulk conductivity can be 
translated, within a Drude picture, into kpT tending to 
a constant value as kp —^ 0. This trend is consistent 



with the experimental observation described above. The 
insensitivity of the conductance to the gate voltage can 
also be expected to occur when puddles are formed with 
a finite (positive or negative) charge density 3^ 37. 3^. 



The analysis presented so far describes the observed 
low temperature gap in transport measurements in terms 
of the features of isolated junctions, neglecting effects 
associated to interference effects between multiple junc- 
tions. We do not expect these effects to change signifi- 
cantly the analysis. Coulomb blockade leads to the sup- 
pression of phase coherence between successive tunneling 
events. The scaling equation ^ is valid for granular ar- 
rays, and transport gaps in these systems have the same 
functional dependence as in single grains [3^ . Hence, in 
a disordered system the transport properties will be dom- 
inated by the junctions with the highest gaps. Note also 
that inelastic cotunneling processes 4^, 41 ] , which influ- 
ence the conductance of single junctions at low temper- 
atures, and which are not considered here, are strongly 
suppressed in junction arrays. 

In conclusion, we find that the gaps observed in con- 
ductance measurements on graphene nanoribbons in Ref. 
[l3 ] can be explained as Coulomb gaps due to the ex- 
istence of internal junctions between graphene islands, 
where the transmission, for all transverse channels, is 
less than one. We identify the gaps observed in trans- 
port measurements as the effective charging energy of 
the islands, renormalized by the charge fluctuations at 
the junctions. We further simplify the model by assum- 
ing that the transport properties can be studied by ana- 
lyzing a single representative junction. The model leads 
to a simple dependence of the gaps on material param- 
eters, such as the carrier concentration, or the width of 
the ribbons. We obtain a reasonable agreement with the 
experiments. 

Charging effects are mostly determined by the geom- 
etry of the system and by the amount of screening of 
the Coulomb interaction. The explanation proposed here 
implies that the gaps observed in transport measure- 
ments should be weakly affected by static disorder, or 
by changes in the electron interference properties, such 
as those induced by an applied magnetic field. On the 
other hand, we expect that charging effects should be 
suppressed by metallic gates at distances from the rib- 
bon which are smaller than the ribbon width. 
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